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Abstract. In this work, applying general results from averaging theory, we find periodic 
^ | 1 orbits for a class of Hamiltonian systems H whose potential models the motion of elliptic 

■ galaxies. 
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1. INTRODUCTION AND STATEMENT OF THE MAIN RESULTS 



Galactic dynamics is a branch of Astrophysics whose development started only around 
sixty years ago, when it was possible to have a view of the physical world beyond the 
integrable and near integrable systems [2] . Even the importance of the analysis of galactic 
potentials, the global dynamics of galaxies is not a simple question and represents a big 
challenge for the researches in the field [2] . Most of the work in the analysis of galaxies is 
^. ■ numerical, in this paper we present an analytical technique, the averaging theory, which 

C"- . allows to find periodic orbits of a differential system. 

In the last years, great quantity of the research on galactic dynamics has been focused 
on models of elliptical galaxies. In most of these models the terms in the potential are of 
even order, so we have adopted this fact in the Hamiltonian system that we are analyzing. 
Another important point that appears in these kind of potentials is that the existence of 
periodic orbits is a useful tool for constructing new and more complicated self consistent 
models. One way to identify periodic orbits is to localize the central fixed points on the 
surfaces of constant energy. In [11J, the authors study the localization of periodic orbits 
and their linear stability for a particular two-component galactic potential. In fact, in our 
days the study of individual orbits in some galactic potentials is a new branch of galactic 
dynamics (see for instance the articles [U EJ |8]) . 

The calculation of particular orbits in some analytical potentials modeling elliptical 
galaxies, indicates that relatively small symmetry breaking corrections can increase dra- 
matically the number of stochastic orbits, showing the importance of the study of pertur- 
bations of simple models |7j . The class of potentials studied in this paper have not chosen 
with the aim of modeling some particular galaxies, our objective is to study systems which 
are generic in their basic properties. 
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In [12], the authors study the galactic potential 

H=^P 2 +P 2 ) + V(X 2 ,Y 2 ). 

These kind of potentials are important in the modeling of elliptic galaxies, as for instance 
we can mention the potentials Vl = log (1 + X 2 + Y 2 /q) and Vc = a/1 + X 2 + Y 2 /q — 1, 
where the parameter q gives the eccentricity of the elliptic galaxy In this paper we deal 
with the Hamiltonian 



H 



(aX 4 + bX 2 Y 2 + cY 4 



(1) 



and its respective Hamilton's equation 
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{4aX 3 + 2bXY 2 ), 
(2bX 2 Y + 4cF 3 ), 

the matrix of the linear part of this system at the origin of coordinates is 

/ 

M = 

V 

with eigenvalues ±z, ±i/q where % = 

In order to obtain periodic orbits for these kind of potentials we will apply averaging 
theory, in this way we re-parametrize the coordinates by the factor y/e for e positive small 
enough (a similar change of coordinates has been used in [9]), that is we do the change 

(X,Y,P X ,P Y ) -> y/i{x,y, Px , Pv ). (3) 

After straightforward computations we get the Hamilton's equations 
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e(4ax 3 + 2bxy 2 ) 
4cy 3 ) 



— — e{2bx 2 y 
Q 



which have the same linear part at the origin than the previous one, the structure of the 
new Hamiltonian is identity with (CQ) in the new variables. Our goal is to study which 
periodic orbits for e = (the unperturbed system) persists for s positive and small enough 
(the perturbed system). 

By the form of the matrix M we observe the necessity to split the analysis for the 
periodic orbits in two cases 
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• q is an irrational number. Here the linear part of system (jlj) has two planes foliated 
by periodic orbits. In the first one the orbits have period 2it, each periodic orbit 
on this plane is of the form 

PO\ = (x cos t + p xo sin t , 0, p xo cos t — x sin t , 0) . 

In the second one, the orbits have period 2wq, each periodic orbit on this plane is 
of the form 

P0 2 = (0,yo cos (t/q)+P yo sin (t/q), 0,p yo cos (t/q) -£/ sin (t/q)). 

• q is a rational number. Here the linear part of system @ has a 4-dimensional 
space filled of periodic orbits of period 2nr if q = r/s with (r, s) = 1, where each 
periodic orbit is of the form 

P0 3 = (a?o cost + p xo sin t, yo cos (st/r) +p yo sin (st/r), 

p xo cost — x smt,p yo cos ((st/r) — ^sin (st/r)). 

When q is an irrational number our main result is: 

Theorem 1.1. For q an irrational number, we have that in every energy level H = h > 
the Hamiltonian system (J2J) has 

(a) at least one periodic solution (X(t),Y(t), Px(t), Py(t) such that when e —> 0, we 
have that (X(0),Y(0), P x (0), P Y (0)) tends to (0,0,0,0); 

(b) at least one periodic solution (X(t),Y(t), Px(t), Py(t)) such that when e — > 0, we 
have that (X(0), F(0), P x (0), Py(0)) tends to (0,0,0,0). 

So, for q irrational, we obtain that in every energy level H = h > the perturbed Hamil- 
tonian system has at least 2 periodic orbits. 

Remark 1.1. We note that the periodic orbits found in the statements of Theorem II .11 are 
in fact degenerate Hopf bifurcations periodic orbits, since they born from the equilibrium 
point localized at the origin of coordinates. Unfortunately we cannot obtain periodic 
solutions when q is a rational number, see Remark 13.11 

The paper is organized as follows. In section 2 we present the theorem from averaging 
theory necessary to prove our main result. In section 3 we give the proof of Theorem II .![ 

2. Some results from averaging theory 

In order to have a self contained paper, in this section we present the basic results from 
the averaging theory that are necessary for proving the main results of this paper. 

2.1. Results from averaging theory. We consider the problem of the bifurcation of 
T-periodic solutions from the differential system 

x'(t) = F (t,x)+eF 1 (t,x), (5) 

where the functions Fq, F\ : K x Q — > M. n are of class C 2 functions, T-periodic in the first 
variable, and Q is an open subset of 1R™. When e = we get the unperturbed system 

x'(t) = Fo(t,x). (6) 
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One of the main assumptions on the above system is that it has a submanifold of periodic 
solutions. A solution of system (jSJ), for e sufficiently small is given using the averaging 
theory. For a general introduction to the averaging theory see the books of Sanders and 
Verhulst [H] , and of Verhulst [15] . 

E 3 ' E 3 

Let x(t, z) be the solution of the unperturbed system (jSJ) such that x(0, z) = z. We 
write the linearization of the unperturbed system along the periodic solution x(i, z) as 

y' = AcF (t,x(t,z))y. (7) 

In what follows we denote by M z (t) some fundamental matrix of the linear differential 
system ([70, and by £ : IR fc x R n_fc — y M. k the projection of R n onto its first k coordinates; 
i.e. £(xi, . . . ,x n ) = (xi, . . .,x k ). 

Theorem 2.1. Let V C R h be open and bounded, and let (3 : C1(V) M n ~ fc be a C 2 
function. We assume that 

(i) Z = {z a = (a,Po(a)) , a e Cl(V^)} C Q and that for each z a e Z the solution 
x(i, z a ) o/ (J6j) is T -periodic; 

(ii) /or eac/i z a G i/iere is a fundamental matrix M Za (t) of (J7J) suc/i i/jai £/ie matrix 
M Zc ^{0) — M~^(T) has in the upper right corner the k x (n — k) zero matrix, and 
in the lower right corner a (n — k) x (n — k) matrix A a with det(A a ) ^ 0. 

We consider the function J 7 : C1(V) — > M fc 

JF(a) = £ (J^ M^{t)F,{t^{t,^)dt \ . (8) 

If there exists a J-'(a) = and det ((dJ r /da) (a)) 7^ 0, t/ien t/iere zs a T -periodic 

solution <p(t,e) of system (j5J) sac/i £/ia£ 99(0, e) — >• z a as £ — > 0. 

Theorem 12.11 goes back to Malkin [TU] and Roseau [13] , for a shorter proof see [3] • 

3. Proof of the main Theorem 

In this section we give the proof of our main result Theorem 11.11 

3.1. Proof of Theorem 11.11 We know that the periodic orbits of a Hamiltonian system 
always appear in cylinders foliated by periodic orbits, each periodic orbit corresponds 
to a different value of the energy h, see for more details pQ. In order to have isolated 
periodic orbits and be able to apply the averaging theory we fix the total energy H = h. 
Computing p x in the energy level H = h we get 

p x = ±^2h -p 2 y /q - x 2 - y 2 /q - e(2ax 4 + 2bx 2 y 2 + 2q/ 4 ). (9) 

The fix value h of the total energy is determined by the initial periodic orbit, which in 
our case for the periodic orbit PO\ it corresponds to h = \{p 2 XQ + Xq), choosing the sign 
+ for p x , and expanding around e = 0we obtain 

rz : z ~ z~~ ax 4 + bx 2 v 2 + cv A . . 
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The equations of motion on the energy level H = (p 2 Q + Xq)/2 are 

ax 4 + bx 2 y 2 + cy 4 



P 2 x -x 2 + x 2 Q - (y 2 + p 2 y )/q 



y = ^, (ii) 
q 

p y = ~ -e(2bx 2 y + 4q/ 3 ). 

In order to apply Theorem 12.11 to system (ITT1) , let 

x = (x,y,p y ), 

F (t,x) = (Jp1 ~x 2 +xl- (y 2 +Py)/q, Py/q, -y/q) 

Fi(t,x) = (_ o^+taVW o _(26x 2 i/ + 4c 2 / 3 



(12) 



The set = {(x, y,p y )\q(pl — x 2 + Xq) — y 2 — p 2 ^ 0} is an open subset of M 3 . Clearly 
the above functions are of class C 2 (Q). The set V of Theorem 12.11 is given by 

V = {z = (xq, 0, 0) : |xo| < p} for some p large enough. 

Let x(i, z) be the solution of the unperturbed system (JBJ) such that x(0, z) = z. The 
variational equations of the unperturbed system along the periodic solution PO\ are 

y , = J D x F (t,x(t,z))y, (13) 

where y is a 3 x 3 matrix. 

The fundamental matrix M(t) of the differential system (fTB"]) such that M(0) is the 
identity matrix of 1R 3 takes the simple form 

(cos t — xq sin (t)/p xo \ 

cos(t/g) sin(t/g) , (14) 

— sin(t/g) cos (t/q) J 

whose inverse is given by 

(Px 1 (Pxo cos t-x sin t) 
cos (t/q) -sin (t/q) 

sin (t/q) cos (t/q) 

An easy computation shows that 

/ 

M _1 (0) - M-\2n) = 2sin 2 (7r/g) sin (27r/g) 

\0 -sin(27r/g) 2sin 2 (7r/g) 

We observe that this matrix has a couple of zeros in the upper right corner of size 1x2; the 
determinant of the 2x2 matrix which appears in the lower right corner is 4 sin 2 (n/q) ^ 
because q is an irrational number. Consequently all the assumptions of Theorem 12.11 
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are satisfied. Therefore we must compute the simple zeroes of the function T defined in 
Theorem 12.11 A straightforward computations shows that 



, a(xn cost + p Xn sin t) 4 
F x t,x t,z = — ^ ^° . / ,0,0 



therefore we get 



^ m(t ,^,. ))= (- ^^ g y ,o,o) 



Let 



f/ \ 1 /" * «Px (^oCost + p X0 sint) 4 2 2 /~7 " 

/i(^o) = 75- / ; : . ^ 2 dt = 3a Pxo(Px + x )/2 = 3ahJ2h - x 2 . 

Zn Jo [Pxo cos t — x sm t) v 

In the last equality we have gotten p xo from the energy relation h = {p 2 X0 + Xq)/2. So the 



solutions of f\ = are xq = ±v2/i, which are simple zeroes. On the other hand we can 
verify that both zeroes generate the same periodic orbit. Then doing the rescaling ([3]) we 
obtain statement (a) of Theorem 11.11 

For the proof of statement (b), as in the previous case we fix the value of the total 
energy as H = {j> 2 , Q + x 2 ] ) /2q = h determined by the initial periodic orbit POi- Computing 
p y from the equation H = h we obtain 

p y = ±^/2qh — qp 2 — qx 2 — y 2 — e(2aqx A + 2bqx 2 y 2 + 2cqy A ), (15) 
we choose the sign + for p y and expand arounf e = getting 

p y = \/2qh — q(p 2 + x 2 ) — y 2 — e— + bx y + y ) ^ 

\/2qh - q(p 2 x + x 2 ) - y 2 

Now, we write the equations of motion on the energy level H = (p 2 o + x^)j2q in the 
order (y,x,p y ), they are given by the system 

\/2qh — q{p 2 + x 2 ) — y 2 ax A + bx 2 y 2 + 2y 4 

y £ 



\/2qh — q(p 2 + x 2 ) — y 2 
x = p x , (17) 
py = -- - e(2bx 2 y + 4q/ 3 ). 

In order to apply Theorem 12.11 to system (TIT}) we are using the same notations and 
definitions (with the obvious changes) than in the previous case. 

Let x(t, z) be the solution of the unperturbed system (jSJ) such that x(0, z) = z. The 
variational equations of the unperturbed system along the periodic solution P0 2 are 

y / = J D x F (t,x(t,z))y, (18) 

where y is a 3 x 3 matrix. 
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The fundamental matrix M(t) of the differential system ( |T8|) such that M(0) is the 
identity matrix of M. 3 takes the simple form 

/ cos (t/ q) - y sin (t/ q) /p yo 
M(t) = cost sint | , (19) 

y — sin t cos t 

whose inverse is given by 

PyJ (Py cos (t/q) - y sin (t/q)) 
M-\t) = I 








An easy computation shows that 





M _1 (0) - M-\2n) = | 2sin 2 (7rg) sin (2rrq) 

— sin (2irq) 2 sin 2 (irq) 

We observe that this matrix has two zeros in the upper right corner of size 1x2; the 
determinant of the 2x2 matrix which appears in the lower right corner is 4 sin 2 (nq) ^ 
because q is an irrational number. Consequently all the assumptions of Theorem 12.11 
are satisfied. Therefore we must compute the simple zeroes of the function J 7 defined in 
Theorem 12.11 

A straightforward computations shows that 

p u W/ ^ ( c (Vo CQS (t/q) + Pvo sin (*/g)) 4 n n 

.Fi(t,x(t,zJJ = — -— — — — — ,U,U 

V Pyo cos {t/q) -y sm(t/q) 

therefore we have 

M (t)F l( t,x(t,z)) - \- {pyoCOs{t/q) _ yoSin{t/qW .0.0) 



Let 



f,. = J_ f 2n f cp yo (y cos (t/q) + p yo sin (t/ g)Y . ^ 
27r io V (Py cos (t/q) -y sm(t/q)) 2 

2. /o*,„ _ .,2 



= 3c^ (p 2 + i/ 2 )/2 = 3chq 2 ^2hq - y 2 . 
In the last equality we have gotten p Xo from the energy relation h = x ° 2q ■ So the 



solutions of /i = are yo = ±^2hq, which are simple zeroes. On the other hand we can 
verify that both zeroes generate the same periodic orbit. Doing the rescaling ([3]) we get 
the statement (b) of Theorem II. 11 

Therefore we have proved that for q an irrational number, in every energy level h > 0, 
the Hamiltonian system (TSJ) has at least 2 periodic orbits, so Theorem 11.11 holds. 
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Remark 3.1. Using the methods of averaging theory studied in this paper, we could not 
obtain any periodic orbit for the Hamiltonian system (F2J) when q is a rational number. We 
have tried to get some information in two different ways, using cartesian coordinates as 
in statement (a) and using a modified kind of polar coordinates in two different planes. 
In the first way we have obtained the variational equations, but unfortunately we could 
not solve them. In the second way we have obtained that one of the equations that we 
must solve for obtain the periodic solutions is identically zero. 
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